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Superlattices (SLs) in monolayer and bilayer graphene lead to extra finite energy and zero energy 
Dirac fermions with tunable anisotropic velocities. We theoretically show that transport in a weak 
perpendicular (orbital) magnetic field allows one to not only probe the number of Dirac points but 
also yields further information about their dispersion. We find that a moderate magnetic field in 
monolayer graphene can lead to a strong reversal of the transport anisotropy imposed by the SL 
potential, due to the SL induced dispersion of the zero energy Landau levels. This effect may find 
useful device applications. For bilayer graphene, we discuss the structure of Landau level wave 
functions and local density of states in the presence of a uniform bias, as well as in the presence 
of a kink in the bias which leads to topologically bound 'edge states'. We consider implications for 
valley filtering and impurity induced breakdown of the quantum Hall effect in bilayer graphene. 



Graphene, a 'Dirac material', exhibits such novel phe- 
nomena as an anomalous integer quantum Hall effect, 
weak anti- localization, and 'zitterbewegung', which stem 
from its low energy excitations being massless chiral 
Dirac fermions.fSS In addition, it holds great potential 
for applications due to its high mobility and the exper- 
imental ability to tune its properties via gating.^^l Early 
experiments aiming to characterize the electronic prop- 
erties of graphene used quantum Hall measurements to 
probe its underlying structure and unequivocally demon- 
strate the relativistic nature of the charge carrier^^li) jj^ 
these studies, the Hall conductivity displayed plateaus 
at atypical values of ±4(e^//i)(A^ -I- 1/2), where N is an 
integer. The origin of these plateaus is well known to 
be a consequence of the electrons in graphene being gov- 
erned by a relativistic Dirac Hamiltonian, together with 
the inherent valley and spin degeneracy. The chiral sym- 
metry of the Dirac Hamiltonian generates a particle-hole 
symmetric Landau level (LL) spectrum - every positive 
energy LL has a conjugate negative energy LL that to- 
gether are responsible for the positive and negative con- 
ductivity plateaus. The Dirac Hamiltonian also supports 
an additional zero energy LL, and leads to the 'half-step' 
offset in the first conductivity plateaus. Finally, the step 
size between each plateau can be attributed to the pres- 
ence of four degenerate Dirac cones labelled by different 
spin/valley indices. Further information about the veloc- 
ity of the Dirac cone, vp, can also be obtained by mea- 
suring the energy gap between the LL, which scales as 
Vf\/~B. Thus, information about LL spectrum and Hall 
conductivity of a system can provide direct evidence for 
the existence of Dirac-like quasiparticles and can be used 
to probe the degeneracy and velocity of the Dirac cones. 

Inspired by studies of superlattices (SLs) in semicon- 
ductor system^, such SLs have been proposed as a route 
to band structure engineering in graphene, paving the 
way to further new physics and applications. Recent 
studies of one dimensional (ID) superlattices (SLs) in 
monolayer graphene (MLG) have demonstrated their re- 



markable ability to not only regulate the Fermi velocity 
but also generate new Dirac cones.'^M^ In seminal work. 
Park et al.,^ showed that a one dimensional (ID) SL in 
MLG could lead to a strong Fermi velocity renormaliza- 
tion in the direction transverse to the SL modulation, 
leading to strong transport anisotropies. Later, it was 
shown that SLs can even lead to nonzero energy Dirac 
point an d multiple zero energy Dirac points near each 
valley.l^Eni These effects stem from the chirality of the low 
energy quasiparticles in MLG. The new Dirac points gen- 
erated by the SL lead to extra zero energy Landau levels 
for SLs in a weak orbital magnetic field. ^ 

For bilayer graphene (BLG), a minimal model of the 
low energy physics has quadratic band touching points 
which can be ga pped out by an electric field perpendic- 
ular to the layers.f^il^ ID SLs in BLG are also predicted 
to display remarkable restructuring of the band disper- 
sion, leading to new Dirac points in the spectrum which 
have anisotropic velocities determined by the SL poten- 
tial, and are perturbatively stable against interactions 
(unlike the quadratic band touching point). For bilayer 
systems, there are two possible types of electrostatic SLs: 
(i) a chemical potential modulation where the potential 
profiles are the same on both layers, or (ii) an interlayer 
bias modulation where the layer potentials are of opposite 
sign leading to a periodically modulated electric field. 

The effect of a chemical potential SL is highly depen- 
dent on the the strength of the modulation — the band 
structure goes through a series of stages as the stren gth of 
the modulation is swept from weak to strong.fi^M] start- 
ing from a weak potential, the quadratic band touch- 
ing point splits into two anisotropic Dirac cones. This 
happens while preserving the net "pseudospin winding 
number" .^23 These Dirac points continually move out 
towards the mini-Brillouin zone (MBZ) befor e becom- 
ing gapped at a certain critical strength.^iSESI Por even 
stronger SLs potentials, the widening of the band gap 
reverts and eventually closes by forming a series of new 
zero-energy Dirac points. The protection of the Dirac 
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points has been shown to be directly relat ed to the chi- 
rality of the low energy BLG quasiparticlesP^'^Sl 

In the case of symmetric ID electric field SLs, four 
anisotropic Dirac cones are generated, two at zero en- 
ergy and two at finite energy. The existence of these 
Dirac points can be see iPSl to aris e from a series of cou- 
pled 'topological' mode^i^l2lllIl2ll that emerge along the 
zero-line where the parity of interlayer bias reverses. In- 
terestingly, these Dirac points are exceptionally robust 
to the size of the bias modulation. In the limit of a long 
period SL, these Dirac points evolve into independent 
topologically confined one-dimensional chiral modes. 

In this paper, we investigate the electronic properties 
of both bilayer and monolayer graphene ID SLs in the 
presence of a uniform magnetic field perpendicular to 
the plane. For the bilayer, we consider two types of 
SLs: a chemical potential modulation and an electric field 
modulation. Similar to early magnetic field studies of 
graphene and graphene SLs'^ ^ careful analysis of the 
Landau levels spectrum and Hall conductivity proffers 
many potential experimental signatures of the underly- 
ing Dirac cones generated by the SL. Special attention is 
also made to electric field modulations where the period 
length is large and decoupled 'topological' modes emerge. 
Although our SL results are most relevant to graphene 
with chiral low energy excitations, the effects explored in 
this paper may also find interesting counterparts in low 
density two dimensional electron gas in semiconductor 
heterostructures, in which honeycomb pattern has been 
successfully imprinted, with a possibility of other kinds 
of SL patterns imprinted in the future.!^ 

We obtain the following key results for ID SLs in MLG 
and BLG in a perpendicular (orbital) magnetic field. (1) 
We show that transport in a weak magnetic field can not 
only be used to probe for the presence of extra Dirac 
points, but may also be used, in clean systems, to ob- 
tain further information about their velocities via the en- 
ergy and degeneracy of higher Landau levels (LLs). (2) 
We show that a moderate magnetic field in monolayer 
graphene leads to a dramatic reversal of the transport 
anisotropy generated by the SL potential alone, an effect 
arising from the SL induced dispersion of the zero en- 
ergy LL. This field tunable transport anisotropy, which 
may find useful applications in monolayer graphene SLs, 
is found to be absent in BLG. (3) Finally, we consider the 
Landau levels in BLG in the presence of a uniform bias, 
and a kink in the bias which leads to ID topologically 
bound states, and the coupling between such modes in 
an array of kinks. We discuss the real space structure 
of the LL wave functions and the local density of states 
which is expected to be relevant to STM measurements 
such as those in Ref. [3^ and EHl and (4) consider possi- 
ble implications for valley filtering and breakdown of the 
quantum Hall effect. 

The organization of this paper is as follows: Section 
H] concerns ID superlattices in monolayer graphene in a 
magnetic field. We begin this section with a brief descrip- 
tion of the band dispersion of ID superlattices in the ab- 



sence of a magnetic field. In the next part, Section YK we 
introduce the low energy effective Hamiltonian and then 
identify the following three magnetic field regimes that 
lead to very different electronic dispersions: (i) weak, (ii) 
intermediate and (iii) high fields. We present and the 
Landau level dispersions of these three regimes in Sec- 
tion [iXl] - |TX3j respectively. We then present the diag- 
onal and Hall conductivity in Section |I B| We find that 
the diagonal conductivity shows strong anisotropy rever- 
sal and Hall conductivity no longer exhibits plateaus, 
when magnetic field is tuned from weak to intermedi- 
ate strength. For realistic SL period, A ~ 300a with 
a — I.42A being the nearest neighbor lattice constant, 
the field strengths corresponding to weak and interme- 
diate magnetic field are about O.IT and 8T, respec- 
tively. Also, for this SL period, the minimum SL poten- 
tial required to generate three Dirac points in one val- 
ley is about 200meV. All these values are experimentally 
achievable in order to see interesting physics. 

Section [H] of this paper examines ID superlattices in 
bilayer graphene and begins with description of the low 
energy theory that describes both chemical potential and 
electric field SLs. We discuss the LL dispersion and trans- 
port for ID chemical potential and electric field SLs in 



Sections II A and II B respectively. 

Then in the Section [ml we consider the real space pic- 
ture of interlayer modulations with dilute kinks and an- 
tikinks that form decoupled soliton modes. We first pro- 
vide an overview of the LL spectrum of uniformly biased 
BLG in Section |III A| and then compare it to the spec- 
trum generated when there is a single isolated antikink in 



the interlayer bias in the beginning of Section IIIB[ Also 
in this section, we present the LDOS and provide key 
signatures that could identify the presence of the kink in 
STM measurements and suggest how the occurrence of 
disorder induced kinks may contribute to the breakdown 
of the quantum Hall effect. Following the discussion of 
a single isolated interlayer bias kink, we discuss a pe- 



riodic array of kink and antikinks in Section III C We 



argue that for £b < A, it is only necessary to consider 
a single kink-antikink pair and observe that regardless 
of the magnetic field strength, the low energy physics is 
governed by the soliton modes which remain robust. Em- 
ploying a simple low energy effective theory to describe 
the soliton modes in a magnetic field, we gain further in- 
sights into the low energy physics of the SL and discuss 
why the zero energy LLs are robust from this perspective. 
Finally, in Section flll D| we consider the coupling between 
the soliton states generated by a kink-antikink pair and 
how this leads to a single ID one-way conducting band 
with definite valley index in each kink separately. 



MONOLAYER GRAPHENE 
SUPERLATTICES 



It is useful to briefly review the effect of the SL 
potentiaP'iSliSIila.t zero magnetic field on the low energy 
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spectrum near K before we present the SL spectrum in 
a nonzero magnetic field. For [7 = 0, where U is the SL 
strength, the spectrum consists of a single Dirac point at 
K, with an isotropic velocity at low energy. The first ef- 
fect of turning on a nonzero U is to make this Dirac cone 
anisotropic, with Vy — Vf, but Vx < Vf- For U — Anvf/X, 
with A as the SL period, one finds = 0. Further in- 
creasing U, with Anvf/X < U < Snvf/X, leads to three 
zero energy Dirac points: one of these continues to be 
located at K, while two new Dirac points emerge which 
are symmetrically split off from the K point along the 
±a;-direction. All three Dirac cones have anisotropic ve- 
locities. It is clearly convenient to measure the strength 
of the SL potential in terms of a dimensionless param- 
eter U, defined via U = {27rv f/X)U, so that tuning U 
with < [/ < 4 leads to the different interesting spectra 
discussed above. A further increase in U leads to a more 
complex spectrum with even more Dirac points. 



A. Landau Level spectrum 

Ignoring spin, the low energy Hamiltonian of pris- 
tine MLG is given by a 2 x 2 matrix at each valley, 
Hq = Vf{sTrx(Tx — T^yO'y), whcrc pseudospin az = ±1 la- 
bels the two sublattices, while the two (decoupled) val- 
leys at ±K = ±47ri:/3 are labelled by s = ±1. Here, 
Vj:=3ta/2 is the isotropic Fermi velocity, with a = 1.42A 
and t = 3eV being the nearest neighbor carbon-carbon 
distance and transfer integral respectively. (We set h=l 
for convenience.) In a uniform perpendicular magnetic 
field, TTj = j — eAj] for B = — and in the Landau 
gauge, the vector potential A — Byx. 

Diagonalizing Hq, we find energies e„ — sgn(n)-y/|n[wcj 
where Wc = V2vp/£b, with £b = l/VeJ3. For s = -1-1 (i.e., 
at valley K), the n^O eigenfunctions are given by 



<f>n,k,+ {x,y) = 



/2L 



V'|n|,fe(y) 



l,k 



(1) 



where L and k are the system length and electron mo- 
mentum deviation from K, both along the a;-direction, 
while for n = 0, 







(2) 



Here, ^n,k{v) is the n-th eigenstate of a (shifted) ID har- 
monic oscillator, 
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centered at j/o = ^^bj and are Hermite polynomials. 
For s = —l (i.e., at — K), the eigenfunctions are given by 
(t>n,k-{x,y) = -i(Jy(t)n,k,+ {x,y). The full low energy LLs 
of MLG are thus (/'„,fc,± (x, ?/)e±^-^-^. 



We now turn to the effect of a periodic ID chem- 
ical potential modulation V{y), with period A ^ a, 
on these Landau levels. The set of eigenfunctions 
(j)n.k,s{x,y)e'^''^^^ , with s — ±1, form a convenient basis 
to study the SL Hamiltonian in a magnetic field. (This 
basis choice is different from the one used by Park, et al, 
and allows us to numerically access a wide range of mag- 
netic fields.) Due to momentum conservation along the 
x-direction, the SL Hamiltonian is diagonal in k. Fur- 
ther, for A ^ a, intervalley scattering is strongly sup- 
pressed. We will therefore assume that the two valleys 
stay completely decoupled. (We focus below on valley K 
with s = -1-1; we expect identical physics around valley 
— K.) With this approximation, the only effect of the SL 
potential is, thus, to induce Landau level mixing. 

To proceed, we need to choose a concrete form for the 
SL potential. For simplicity, we set V{y) — ^ cos (^), 
although our results can be easily generalized to other 
(e.g., step- like) SL potentials by including multiple 
Fourier components. For numerical computations, the 
following integral identitjJ^ for harmonic oscillator states 
proves useful: 



dy cos J ^n,k{y)'^n+7n,k{y) 



(n + m) ! 



1/2/ 



V2tt£i 



exp 



A / V 

' 2nyo rnw" 

X : 



(4) 



where m > 0, and Ln"'' is the generalized Laguerre poly- 
nomial. Using this identity, we simplify and numerically 
compute the matrix elements of the full Hamiltonian, re- 
taining up to 3000 Landau levels, and diagonalize this to 
obtain the spectrum of the ID SL in a magnetic field. 

In order to study the effect of the magnetic field on 
the ID SL in graphene, with U ^ 0{1), it is useful to 
consider three regimes for the magnetic field. 

(i) Weak field: This regime corresponds to having 
hujc ^ U, where the Landau level spacing is much smaller 
than the SL amplitude, so that 2Tr£B/X ^ 1. In this 
regime, the magnetic field may be viewed as effectively 
'probing' the zero field SL excitations. 

(ii) Intermediate field: In this regime, huc ~ U, 
which means 2t:£b/X ~ I, so that the SL potential and 
the magnetic field have to be treated on equal footing. 

(iii) Strong field: Here, tiu>c ^ U or, equivalently, 
2it£b/X <^ I. In this regime, the SL potential only 
weakly perturbs the Landau levels of pristine graphene. 

Fig. [T] shows the spectrum of the graphene SL in differ- 
ent field regimes for SL strengths U—2TTVf/X (or U — 1) 
and QiTVf/X (or [7 = 3). This allows us to contrast 
the behaviour of the spectrum of the SL in a magnetic 
field without or with extra Dirac points being present at 
zero field, and to explore consequences for quantum Hall 
physics and transport. 
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(a) 



(c) 



(b) 




tred at K has a slightly different average velocity ^v^Vy 
compared with the two cones which are symmetrically 
split off from K along ±i. This degeneracy breaking re- 
sults in the Landau levels at nonzero energy becoming 
weakly split, as is most clearly seen for the first two ex- 
cited Landau levels (at positive or negative energy, i.e., 
with n = ±1,±2). We have numerically determined 
and Vy for each of the three Dirac points and found good 
agreement between the energy levels obtained on this ba- 
sis of having Dirac fermions with two different average 
velocities, and that obtained directly numerically. 

At higher energies, iJ/wc 2 for [/ = 1 or iJ/wc > 1 
for U = 3, the spectrum begins to deviate from this sim- 
ple behavior expected for a linear Dirac spectrum. This 
deviation results from curvature in the dispersion, which 
appears upon going beyond the linearized approximation. 



FIG. 1. (Color online) Landau levels of monolayer graphene 
SL for different (dimensionless) SL strengths U, and magnetic 
fields B. The spectrum is shown for weak field {£b ~ 2A, top 
panels) and intermediate field {£b = 0.2A, bottom panels). 
Left panels (a,b) correspond to C/ = 1 which supports a single 
anisotropic zero energy massless Dirac fermion. Right panels 
(c,d) correspond to U = 3 which supports three zero energy 
massless Dirac fermions with anisotropic velocities - the weak 
field zero energy LL thus has three times as many states for 
f/ = 3 as it does for [7 = 1, while the n = ±1, ±2 levels have 
degeneracy splitting in weak field due to the Dirac fermions 
having two different mean velocities. For £b <^ X (not shown), 
the LLs closely resemble that of pristine graphene. See text 
for a detailed discussion of the Landau level structure. 



1. Weak field regime 

When the magnetic field is weak, £b = 2A (top panels 
in Fig. [T]), we find that the energy spectrum barely de- 
pends on the value of fc, or equivalently, ya. This is due 
to the fact that when magnetic length £b is larger than 
the SL period A, the matrix elements of the Hamiltonian 
do not depend on the center of the LL wave functions, 
which yields flat bands. Equivalently, in this regime, the 
magnetic field may be viewed as effectively 'probing' the 
structure of the zero field SL dispersion leading to Lan- 
dau levels which depend on the nature of the Dirac spec- 
trum at low energy. 

For U — 1, the low energy spectrum of the SL con- 
tains a single anisotropic Dirac point at zero energy. 
For an anisotropic Dirac cone described by an effective 
Hamiltonian H — v^kxix + ^y^yCy, the LLs are given 
by Sn = sgn{n) y^2\n\vxVy / £ B ■ Since the SL renormal- 
izes Vx < Vf, but leaves Vy = Vf, the Landau levels at 
weak field resemble those of pristine graphene, but with 
a renormalized lower effective velocity ^VxVy < Vf. 

For U = 2>, the low energy spectrum of the SL contains 
three anisotropic Dirac points at zero energy, so that the 
zero energy Landau level has three times the degeneracy 
of the case with [7 = 1. Further, the Dirac cone cen- 



2. Intermediate field regime 

At intermediate fields, for £b = 0.2A, the spectrum 
at low energy is most simply understood as arising from 
the SL potential inducing a strong dispersion to the Lan- 
dau levels. In simple terms, if we assume that the state 
labelled by momentum k, or equivalently position j/qj 
have an energy which is modulated by the SL poten- 
tial, we expect a periodic modulation of this energy with 
period A and amplitude proportional to the SL ampli- 
tude U . The behaviour of the low energy Landau levels, 
n = 0,±1,±2, as seen from the lower panels in Figjl] 
is consistent with this scenario, with the modulation fol- 
lowing the cos(27ry/A) form of the SL potential and the 
modulation for = 3 being roughly thrice as strong as 
the modulation for U = \. We are not able to under- 
stand the dispersion of the higher energy Landau levels 
in such simple terms, although our numerical computa- 
tion reveals that these are also strongly dispersing. 



3. High field regime 

For very strong magnetic field, the Landau level struc- 
ture of pristine graphene is recovered. Here, only one 
zero energy level exists at the Dirac point, and other en- 
ergy levels follow the square root relation. This is simply 
because in such a strong magnetic field, the SL is just a 
perturbation and can only give rise to a small modulation 
of the LLs following our argument at intermediate field. 
From a perturbative point of view, the energy corrections 
up to first order to the LL energies are given by 

Ai?(i'= J dy4,Uk,y)V{y)Mk,y), (5) 

which, upon using the identity Q, gives a sinusoidal 
dependence on the center position of LL wavefunctions. 
Thus, even in a strong magnetic field, the energy spec- 
trum is not dispersionless but has a spatial modulation 
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FIG. 2. (Color online) Diagonal dc conductivities of mono- 
layer graphene SL for different (dimensionless) SL strengths 
U, and magnetic fields B. The conductivity is shown for weak 
field {£b ~ 2A, top panels) and intermediate field {£b ~ 0.2A, 
bottom panels). Left panels (a,b) correspond to f/ = 1, and 
right panels (c,d) correspond to U — 3. The conductivities 
show strong anisotropy when magnetic field strength is tuned 
- for weak field (a,c), ayy is larger than a.^^, which is a conse- 
quence of the Fermi velocity renormalization in the absence of 
magnetic field; for moderate field (b,d), the anisotropy is re- 
versed, since Vx acquires intra-LL contributions, as explained 
in the text. For ^ A (not shown), result for pristine 
graphene is recovered and the transport is isotropic in both 
directions. 



following the SL; however the ratio of the amplitude of 
this modulation to the Landau level spacing, U/uJc, is 
extremely small in the high field regime. 



B. Diagonal and Hall Conductivity in MLG SLs 

Once we have the eigenvalues and eigenfunctions for 
the superlattice in a perpendicular magnetic field, both 
ac and dc conductivities can be calculated directly by 
Kubo formula,!^ 



crjj(w) = 



1 



{ak\vi\/3k) {/3k\vj\ak) 

En — Efj — LU — iT 



f{E^) - f[Ep) 

Ea — Ejj 



(6) 



Here, we have set F = 10"'^ x 27ru//A as the Landau level 
broadening, Ea{yo) and \ak) are the a-th eigenvalue and 
the corresponding eigenstate of the system which can be 
expanded in the basis of \nk), where 4>n{k,y) = {y\nk) 
is the LL wavefunctions for pristine graphene. Vi is the 
velocity operator in i-direction and Vi = vpai, where ai 
is the Pauli matrix. Using Eq. ([T]), ([2]), we can get the 
expression for the matrix elements of velocity operators 



FIG. 3. (Color online) The dc Hall conductivity of mono- 
layer graphene SL for different (dimensionless) SL strengths 
U, and magnetic fields B. The conductivity is shown for weak 
field {£b ~ 2A, top panels) and intermediate field {£b = 0.2A, 
bottom panels). Left panels (a,b) correspond to f/ = 1, and 
right panels (c,d) correspond to f/ = 3. For weak field (a,c), 
the Hall conductivity shows well-defined plateaus, as a con- 
sequence of nearly flat energy bands. For intermediate field 
(b,d), the energy bands become dispersive and the Hall con- 
ductivity no longer shows step-like structure. However, for 
weak SL (b), the energy bands are not fully overlapped. Hall 
conductivity still shows small plateaus when chemical poten- 
tial falls between two bands, and the value of a^y changes by 
one between adjacent steps, as expected from Dirac physics. 
For <C A (not shown), result for pristine graphene is recov- 
ered and Hall conductivity is constant between adjacent LLs 
and changes by one when chemical potential crosses an LL. 



from 



and 



{mk\vx\nk) = VF[sgn{m)S\m\-i,\n\ 
-Hsgn(n)(5|„|j„|_i], 



(7) 



{mk\vy\nk) = [sgn(TO)(S|„|_ij„| 

-sgn(n)(5|„|j„l_i], (8) 

by expanding \ak) into \nk). Note that {ak\vy\ak) — 
is always true for any state. 

Results for dc diagonal conductivities as function of 
chemical potential /z are shown in Fig. [2j This can be 
done by setting the frequency uj to zero in Eq. ([6]), and 
only the real part of the conductivity tensor is nonzero. 
In weak magnetic field, the conductivities show strong 
anisotropy, with ayy larger than a^x, which is a conse- 
quence of the Fermi velocity renormalization in the ab- 
sence of magnetic field (see Fig. [2] (a) and (c)). Since 
{ak\vy\ak) = and {ak\vx\cek) ~ because of the fiat 
band structure, the major contribution to the diagonal 
conductivities comes from off-diagonal matrix elements, 
{ak\vi\f3k) with a ^ (5. Numerically, we have observed 



6 



that matrix elements of Vy is always larger than those 
of Vx, which gives rise to the anisotropy in the weak 
field. In intermediate magnetic field, conductivities still 
show anisotropy, but with a^x significantly larger than 
tJyy (see Fig. |2](b) and (d)). This is because Vx has ac- 
quired diagonal matrix element, {ak\vx\ak) = dEa{yo = 
k£g)/dk 7^ since the energy spectrum is dispersive, 
while Vy still lacks this contribution. Notice the posi- 
tions of the conductivity peaks of ayy exactly correspond 
to the minimum and maximum of the energy band, where 
the density of states is the largest. For axx, however, the 
conductivity is minimum at the band edge, since the av- 
erage of the velocity operator, (vx), is zero. Therefore, 
the intra-LL contribution to axx is the smallest at the 
band edge. For weak SL potential, drop to zero 

when there is no overlapping LLs, while in a strong SL, 
axx always show dispersive transport property. In strong 
magnetic field (not shown here), where the Landau levels 
structure of pristine graphene is recovered, the conduc- 
tivities become isotropic (see, for example, Ref. [33]). In 
this case, the SL is merely a perturbation to the magnetic 
field and thus should have minor effect on determining 
the magnetotransport properties. 

The dc Hall conductivity is shown in Fig. [Sj For weak 
magnetic field (Fig. |3](a) and (c)), the Hall conductivity 
shows well-defined plateaus, as a consequence of nearly 
fiat energy bands. The values of Hall conductivity around 
Dirac points are ±l/2(e^//i) in weak SL {U = 1) and 
±3/2(e^//i) in strong SL {U — 3). This result resembles 
the anomalous half integer quantum Hall effect in pris- 
tine graphene and the Hall conductivity triples due to 
the existence of three Dirac points in a strong SL. Mov- 
ing away from the Dirac point, we can observe quantum 
Hall plateaus with higher conductivities, and the value 
increases by 1 each time the chemical potential crosses an 
LL. For intermediate ma gnetic field (Fig. [3](b) and (d)), 
there is no longer well defined plateaus due to the dis- 
persive energy spectrum. However, for weak SL, the LLs 
are not overlapped with each other. If chemical potential 
falls between two LLs, a small plateau can still show up, 
with the value expected from Dirac physics. When mag- 
netic field becomes strong enough as the LL structure for 
pristine graphene is restored. Hall conductivity will show 
anomalous half integer quantum Hall plateaus. 

Fig. |4] shows the ac conductivities of graphene SLs 
in an intermediate magnetic field. For weak and strong 
magnetic fields, the results resemble those of pristine 
graphene,'^ since in both cases the LLs are nearly fiat 
and the real part of the conductivities show strong peaks 
when photon frequencies exactly correspond to the en- 
ergy differences between two LLs. In an intermediate 
magnetic field, the result is complicated by the disper- 
sion of LLs. At low frequencies, there can be optical 
transitions in a range of photon energies, and the real 
part of diagonal conductivities is maximum at the band 
edge where the DOS is also maximum. At high frequen- 
cies, the LLs become less dispersive and peaks will show 
up. These results can be linked with graphene's unusual 
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FIG. 4. (Color online) The ac conductivity of monolayer 
graphene SL for different (dimensionless) SL strengths U , in 
intermediate magnetic field, (.b ~ 0.2A, with fi = 0.2ajc. Left 
panels (a,b,c) correspond to C/ = 1, and right panels (d,e,f) 
correspond to (7 = 3. 



magnet o-opti cal properties, for example, giant Faraday 
rotation.ESEH While the anisotropy in the diagonal con- 
ductivities can lead to anisotropic rotation angles for inci- 
dent waves with different polarization plane, this effect is 
actually quite small and hard to observe experimentally. 



II. BILAYER GRAPHENE SUPERLATTICES 

From the previous section, we can see that the phe- 
nomenon resulting from the existence of Dirac points is 
most manifest when magnetic field is weak compared 
to the SL strength. A very special feature of bilayer 
graphene SL is that anisotropic Dirac points can be gen- 
erated in the energy spectrum, where originally only 
quadratic band touching points is present. Our numer- 
ical results also indicate that in intermediate magnetic 
field, the LLs are only slightly dispersive, compared to 
monolayer graphene. Therefore, in this section, we will 
focus on the properties of bilayer graphene SLs in weak 
magnetic fields using a two-band effective Hamiltonian. 
We hope that our results can provide means to the deter- 
mination and characterization of these anisotropic Dirac 
points. 

The low energy physics of bilayer graphene near the 
Brillouin zone corners at ±K is, in the presence of a 
perpendicular magnetic field, described by the following 
two-band effective Hamiltonian, 



Hn = — 











(9) 



with s = ±1 corresponding to ±K as in the case of 
monolayer graphene. Similar to the single layer case, we 



FIG. 5. (Color online) Energy spectrum of BLG subject to 
a chemical potential SL and a weak perpendicular magnetic 
field {£b = 2A). For both panelss, A = 100a, where a = 1.42A. 
In the left panel, Vb = O.Olt, two anisotropic Dirac points 
are generated in the spectrum in the absence of magnetic 
field, and the zero energy LLs are doubly degenerate. In the 
right panel, Vb = 0.04t, no Dirac points are present, and no 
zero energy LL exists. In both panels, the energy bands are 
nearly flat. For both intermediate flelds (not shown here), the 
energy bands are all slightly dispersive, in contrast to MLG 
SLs where energy bands overlap. This can be seen from the 
diagonal conductivities in Fig. [6] 



have replaced the momentum operator with its canonical 
counterpart to take into account of the magnetic field ef- 
fect. In the following, we will work with the same gauge 
choice, A = Byx. For s = +1, the eigenvalues and the 
corresponding eigenvectors of the above Hamiltonian are 



sg 



</>n,fc,+ (a;,?/) = 



■n{n)^\n\{\n\~l)u:l/t^ 



V'|n|,fe(y) 

-sgn(n)V;|„|_2,A;(2/) 



(10) 



with |n| > 2. In addition, there are two zero energy 
solutions. 



0, 

£q = 0, 



7l 



V'i,fc(y) 






(11) 



For s = —1, the corresponding eigenvectors are given 
by <l>n,k-{x,y) = (Jx(t)n,k,+ {x,y). The fuU low energy 
LL wavefunctions thus take the form 0n.fc,±e**^"=^. and 
these serve as a good basis to study the magnetic field 
effect of bilayer graphene SLs. 

The SL can be modelled by assuming spatially varying 
potential profiles on each layer, 



VliK) 
1^2 (x) 



(12) 



where different combination of Vi and V2 gives rise to 
different types of SL, e.g., V^i(x) = V2(x) for chemical 
potential SL and Vi(x) = — V2(x) for electric field SL. In 
following, we will only consider ID SLs, where for proper 
choice of SL strength, various anisotropic Dirac points 
can be generated in the spectrum. 



J 2 



FIG. 6. (Color online) Diagonal dc conductivities of bilayer 
graphene SL for different SL strengths Vb, and magnetic fields 
B. The conductivity is shown for weak field {(.b ~ 2A, top 
panels) and intermediate field {(.b ~ 0.2A, bottom panels). 
Left panels (a,b) correspond to Vb = O.Olt, and right pan- 
els (c,d) correspond to Vo = 0.04t. The conductivities show 
anisotropy, where a^x is always larger than Uyy, in contrast 
to anisotropy reversal in MLG SLs. 



A. Chemical potential superlattice 

For ID chemical potential SLs, we will take V^i(x) = 
V2(x) — V{y) — ^ cos(27ry/A). In previous work,'^^ 
we have shown that for ID chemical potential SL, two 
anisotropic Dirac points will derive from the original 
quadratic band touching points as long as the SL strength 
is smaller than a critical value. This results from the 
chiral symmetry protected level crossing in the SL di- 
rection and level anti-crossing in other directions. Once 
SL strength increases beyond this critical value, Bragg 
reflection will open up a gap in the energy spectrum on 
the mini Brillouin zone boundary. This behavior has also 
been confirmed by Tan et aZ.^ 

Results for energy spectrum of chemical potential SL 
in a weak magnetic field are shown in Fig. [5j In the 
left panel, we have chosen Vq = O.Oli, where t = 3 eV 
is the intralayer nearest neighbor transfer integral. For 
this value of SL potential, our previous calculation has 
shown that there will be two Dirac points generated in 
the mini Brillouin zone. Here, similar to the single layer 
case, the energy spectrum shows flat band structure. At 
the Dirac point, there are two nearly degenerate zero 
energy levels. Right above and below the Dirac point, 
we can also observe two sets of doubly degenerate energy 
levels. These are the n — ±1 LLs of each anisotropic 
Dirac cone. Moving further away from the Dirac point, 
no such degenerate levels are present. This is because 
these Dirac points only have linear dispersion at rather 
low energy. In the right panel, Vq = 0.04t, energy gap 
opens up and there exist nonzero energy levels near the 
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FIG. 7. Evolution of low lying energy levels in a chemical 
potential SL as a function of SL potential strength U , with 
Ib = 2A, 2/0 = 0, A = 100a where a = 1.42A. 



Dirac point. Energy bands in intermediate field regime 
(not shown here) are only slightly dispersive, where the 
amplitudes of dispersion are extremely small compared 
to the LL spacings. 

Fig. [6] shows the dc diagonal conductivity of chem- 
ical potential SLs. Similar to the single layer case, we 
can observe the anisotropy in the conductivity. However, 
there is no anisotropy reversal as magnetic field strength 
is tuned. In weak field, the transport anisotropy is de- 
termined by the anisotropic Dirac cones. From Ref. 1151 
Vx — 2vy for emergent Dirac cones. Therefore, the con- 
ductivity in the x direction, axx, should be larger than 
ayy in a weak magnetic field. For intermediate magnetic 
field, due to the dispersion of the energy bands, the aver- 
age velocity in the x direction is not zero, (vx) 7^ 0. On 
the other hand, (vy) is always equal to zero. This means 
that axx will acquire intra-LL contributions, while ayy is 
mainly determined by inter-LL contributions and is thus 
small compared to (Txx- 

We can further demonstrate how the flat energy lev- 
els evolve as the strength of the SL potential varies, as 
shown in Fig. [7] Here, we have fixed the magnetic field 
strength by setting = 2A and also fix the center posi- 
tion of the wavefunctions to be at 0. As the SL poten- 
tial tends to zero, where the problem reduces to pristine 
bilayer graphene in a perpendicular magnetic field, we 
can see the energies of LLs will follow the pattern of Eq. 
( 10 1. As SL potential increases, the physics gradually be- 



comes dominated by anisotropic Dirac cones, which can 
be seen from the appearance of doubly degenerate levels 
at nonzero energies. This crossover from nonrelativistic 
to relativistic behavior can be qualitatively understood 
by considering the competition between the characteris- 
tic energy scales in these two regimes, fn pristine bi- 
layer graphene, the low energy excitation are massive 
electron with an effective mass m* = t±/2vp. In a 



yn<^> 
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FIG. 8. Energy spectrum and dc diagonal conductivities of 
electric field superlattice, in a weak magnetic field, with (.b = 
2A, U = 0.04t, A = 100a where a = 1.42A. Results for 
intermediate magnetic fields are similar to those of chemical 
potential SLs and are not shown here. 



magnetic field B, the characteristic energy scale is the 
cyclotron frequency, = eB/m*c = l/m*t\. On the 
other hand, the anisotropic Dirac points generated by SL 
have anisotropic Fermi velocities Vy = \/2A|t/(Q)|/7r and 
Vx = '2vy, where Q = y2Ti/\}^ Therefore, the charac- 
teristic energy scale associated with these Dirac points is 
= ^IVxVyl is- We expect to see a smooth crossover 
as these two energy scales are comparable to each other. 
A rough estimate shows that the crossover should occur 
around V ~ 0.002i, which is quite close to the value that 
we have observed. 

Further increasing the SL potential, the doubly degen- 
erate zero energy levels become gapped and all levels are 
pushed away from Dirac point. Surprisingly, at rather 
strong SL potential, U ~ 0.22t, zero energy LLs appear 
again, and all of the higher energy levels become dou- 
bly degenerate. This phenomenon can be understood 
from the result of Tan et al'^ As it has been shown, 
for a chemical potential SL, when SL potential is strong 
enough, anisotropic Dirac cones will show up again in 
the energy spectrum, which naturally leads to the zero 
energy LL in the presence of a magnetic field. Accord- 
ing to Ref. iilG;, for certain values of SL potential, there 
can be four Dirac points in the spectrum, which is not 
present in our calculation. Of course, for strong SL po- 
tential, the two-band description of bilayer graphene is 
no longer valid. However, we have reproduced this result 
in a four-band effective Hamiltonian approach. In strong 
chemical potential SL, there can exist four-fold degener- 
ate zero-energy LLs, and the degeneracy reduces to two 
upon further increasing the SL strength. 



B. Electric field superlattice 

For ID electric field SLs, we will take T4(x) = 
— V2(x) ~ V{y) — ^ cos(27r?//A). Following the same 
procedure as in the previous section, we have calculated 
the energy spectrum and the dc conductivities for an elec- 
tric field SL in a magnetic field. 

In a ID symmetric electric field SL, there are always 
two zero energy Dirac points present in the spectrum. 
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FIG. 9. Evolution of low lying energy levels in an electric field 
SL as a function of SL potential strength U, with £b ~ 2X, 
yo = 0, A = 100a where a = 1.42A. 



which results from the coupUng of ID zero modes at 
kink/antikink of the SL potential profile. Earher, we have 
developed an effective "wire" model to describe these 
Dirac points and demonstrated that they are topological 
stable if the "inversion" symmetry V{y + A/2) = —V{y) 
is preserved. This implies when a weak magnetic field 
is applied, doubly degenerate zero energy levels should 
appear at the Dirac point. Indeed, as it can be seen from 
Fig. |8] and [9j these two levels always stay at zero energy, 
independent of the SL potential. Moreover, for the SL 
potential we choose, Vq — 0.03t, nearly degenerate levels 
appear even at nonzero energies (Fig. [s]). These corre- 
spond to the LLs derived from anisotropic Dirac points, 
up to n = ±4. From Fig. [9j it is more clear that at strong 
SL potential, physics is strongly dominated by the Dirac 
points, where higher energy levels become doubly degen- 
erate and resemble the higher LLs of the Dirac cones. 
When SL potential is weak, equally spaced LLs are re- 
covered, as in the chemical potential SL, which also indi- 
cates a nonrelativistic to relativistic crossover at certain 
SL strength. Remarkably, and different from the chemi- 
cal potential SL case, the relativstic behavior survives to 
higher energies as SL potential increases, which means 
the linear approximation description of Dirac cones works 
in a larger energy range. This is consistent with earlier 
result .f^S Consider a potential profile with a kink of bias 
reversal. When the bias is large, the conduction and va- 
lence bands are widely separated and a pair of linearly 
dispersed zero modes will traverse the gap. Now, when 
these kinks and antikinks are arranged periodically, Dirac 
cones will result from the coupling of these modes, and 
the Fermi velocity along SL direction inherits its value 
from the freestanding zero mode.'^ Therefore, as the SL 
potential increases, the energy range where the Dirac 
cone approximation is valid also increases, which leads 
to the robust relativistic physics at large SL potential. 



III. REAL SPACE PICTURE FOR 
INTERLAYER BIAS MODULATIONS IN 
BILAYER GRAPHENE 

As suggested in an early study by the present 
authoral^, it is extremely useful to view an interlayer 
bias modulation as establishing a series of 'kink' and 'an- 
tikink' modes along the zero-lines where the interlayer 
bias reverses polarity. Within the effective 2x2 low- 
energy description of BLG, the interlayer bias can be 
regarded as a mass generating term and the zero-lines 
represent boundaries where the sign of the mass changes. 
As consequence, along a single isolated kink (antikink), 
a pair of two right (left) moving 'topological' modes at 
K and corres ponding tw o left (right) moving modes at 
the other —J^SMMMM^ The typical length scale over 
which these kink-antikink modes are confined, I, varies 
from - lOOd for AV = O.Oli to < 50 lattice sites for 
AV — O.lt, where d — ^/3a is the length the primitive 
lattice vector. For the remainder of this section we set 
d=l. 

Accordingly, an interlayer bias modulation can be 
thought of as a periodic array of coupled kink-antikink 
solitons. With this viewpoint, and using the symme- 
try properties of the kink and antikink modes, Killi et 
al!^ derived an effective low-energy hamiltonian describ- 
ing how anisotropic Dirac cones precipitate precisely at 
the band crossing point between the kink and antikink 
modes when the soliton modes overlap and couple. With 
this same perspective, further intuition into the quantum 
Hall effect in the SL system can be obtained by studying 
the properties of ID kink states in a magnetic field. 

Aside from its relevance to SLs, a study into the prop- 
erties of the kink/antikink states also provides valuable 
insights into disorder BLG samples in the quantum Hall 
regime. Various sources of disorder are known to generate 
a random electrostatic potential landscape. Throughout 
these samples, ID kink-states are expected to percolate 
along precipitous fluctuations that reverse the parity of 
interlayer bias. If the percolation networks are well ex- 
tended, these kink states would contribute to the break- 
down of the quantum Hall effect^^. Recent STM mea- 
surements of BLG samples in the quantum Hall regime 
also indicate that even when the sample is uniformly bi- 
ased by external gates, the disorder is still strong eriough 
to locally reverse the polarity of the interlayer biaJ^SMl, 
These recent measurements provide an additional moti- 
vation to study the properties of kink states in a mag- 
netic field and, specifically, to examine how they modify 
the local density of states (LDOS). 

In this section, we consider interlayer bias modulations 
in quantum Hall regime from the perspective of the kink- 
antikink modes. We study the dilute limit where the kink 
and anti-kink states are well separated and confined to 
the zero- lines (i.e. A 3> 0- Within this regime, the low 
energy states of the system are completely dominated by 
midgap soliton modes. We start by first reviewing the 
effects of a magnetic field on uniformly biased BLG in 
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section IIII A| and then non- uniformly biased BLG with a 



(a) 0.1 



single isolated anti-kink in the potential profile in section 
|IIIB| . In addition to studying the dispersion relation, we 
examine how the LDOS is modified in the proximately of 
a kink. Next, we consider an array of decoupled kink and 
anti-kink modes in Section [lII C| . For magnetic fields such 
that is < A, it is only necessary to understand the dis- 
persion relation of an isolated kink/anti-kink pair as this 
is sufficient to describe the entire band structure. This 
allows us to consider a zigzag ribbon of width Nceii with 
one period of modulation (i.e. a pair of kink anti-kink 
states) and implement a Peierl's substitution in the full 
tight-binding Hamiltonian. For weaker magnetic fields 
such that £b ?J A, it becomes necessary to include a 
larger number of kink/anti-kink states to describe the 
high energy features of this system. Despite this, fur- 
ther insights into the low energy modes can be made 
by combining the knowledge obtained from the study a 
single kink/anti-kink pair and by using a simple low en- 
ergy effective theory. Aside from this, in section HID we 



demonstrate that a coupled pair of kink/anti-kink states 
opens asymmetric bandgaps in each valley that can be 
precisely controlled using a magnetic field. This effect 
may provide a viable route to developing a switchable 
valley filter. 
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FIG. 10. Landau Level spectrum of graphene strip in the pres- 
ence of (left column) a uniformly interlayer bias and (right 
column) with an anti-kink bias profile with AV = O.OSt, 
N^M = 300. (c,d) n=0 and n=l LLs. Solid red (blue) 
lines are bulk states localized in the upper (lower) half of 
the sample and red (blue) dotted lines are the upper (lower) 
edge states. Black dotted lines connecting the bulk LLs are 
anti-kink states. We have chosen a strong magnetic field of 
B = 120r in order to make the relevant features easily dis- 
cernible and note that the spectrum is qualitatively identical 
at lower magnetic fields. 



A. Uniform Bias 

Although it has been discussed extensively in the liter- 
ature, we review and emphasize a few important points 
about the LL structure in uniformly biased BLG. Details 
about the LL structure in biased BLG can be found in 
Ref. mi [221 1551 - H51 here we only briefly summarize some 
of the features as they pertain to our current discussion. 

In the Landau gauge that preserves translation sym- 
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metry along the zigzag direction, inspection of Fig 
(a) shows that there are two distinct energy level spec- 
trums in the dispersion relation, one about each valley. 
Each energy level consists of two well-defined regions: a 
flat, non-dispersive region consisting of bulk LLs and dis- 
persive edge stateJSIS] (note, the edge state dispersion 
for the n = and n = 1 are unique in that they con- 
tained non-dispersive regions, as seen in Fig. 10 (a,c)). 
The wavefunctions of the non-dispersive LL states are 
well confined within the bulk and, starting from states 
on the far left side of the level and ending on the right 
side, move from one edge of the sample to the other. 

A comparison of the spectrum about the K-points 
clearly shows that the interlay er bia s breaks the valley 
degeneracy of the LL spectrumJ^SBol^ Although the val- 
ley degeneracy is broken, the states in each valley are 
connected by a symmetry that relates the low energy 
Hamiltonian of opposite K-point. If we define the op- 
erator V to correspond to the interchanging oi A ~ B 
sublattices combined with layer exchange (i.e. Ai o B2 
and Bi o A2), it is a simple matter to show that 



r^HK{-AV)r = i?K'(-Ay). Thus, upon interchang- 
ing the sites labelling, the quasiparticles in each valley are 
governed by the same low energy Hamiltonian but with 
the effective relative parity of the interlayer bias reversed, 
leading to the observed valley degeneracy breaking. 



Inspection of the band structure shown in Fig. 10 (c) 
reveals that in the presence of an interlayer bias, the zero 
energy n = and n = 1 LLs are shifted to finite energy 
and their degeneracy is lifted. Note, the energy shift of 
these two these states is positive in one valley an nega- 
tive in the other. A second crucial observation about the 
n = and n = 1 Landau states (and to a lesser extent 
the higher energy states) is that the wavefunctions are 
strongly localized to one of the two layers determined by 
its valley index. Such layer polarization is directly ob- 
served in STM measurements and information about the 
LL spectrum, valley index and local interlayer bias can be 
obtainecf^. Below, we discuss further how the presence 
of a kink and anti-kink affects the LDOS and give rise to 
clear signatures observable in STM measurements. 



B. Single Kink or Antikink in the Bias 

In the absence of a magnetic field, an isolated kink in 
the interlayer bias generates two unidirectional dispersing 
subgap bands in one valley and two oppositely dispersing 
bands in the other valley related by time reversal. An 
anti-kink generates similar low energy bands, although 
the velocity of these modes is reversed in each valley (see 
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Fig. [To] and Ref. [55] for details). 



When a large magnetic field is introduced, the energy 
spectrum resembles that of the LL structure of uniformly 
bias BLG, but with further distinctive features. Consider 
the energy spectrum about a single valley. Comparison 
of Fig. 10 (c) and (d) shows that instead of each en- 
ergy level having a single non-dispersive bulk LL, each 
LL breaks into two flat regions shifted in energy, which 
are connected through a dispersive mode. States in these 
two flat regions correspond trivially to the bulk Landau 
states that would have been generated by having a uni- 
form interlayer bias of the same parity as the local bias. 
Hence, these two regions are composed entirely of Landau 
states that are well localized in the bulk and respond only 
the local interlayer bias. As in the previous case, far way 
from the K-point the bulk states evolve into edge states, 
and are again sensitive only to the local interlayer bias 
at that edge. 

However, it is for the n = and n ~ 1 LLs that the 
presence of the kink states are directly observed. Starting 
from a bulk state in the n = (n = 1) Landau level on 
one side of the interface and tracking the states to the 
other side of interface (moving left to right in momentum 
space), we observe the following sequence. The bulk n — 
(n = 1) LL states begin to continuously evolve into 
dispersive kink states as the guiding center approaches 
the zero-line and then emerge from the interface on the 
other side as bulk states in the n — I {n — 0) LL. 

Although other dispersive states connecting the high 
energy bulk LL on either side of the interface are also 
present, they are not of topological origin. These states 
are simply LL states that becomes dispersive as their 
wavefunction begins to overlap with the region of spa- 
tially varying potential. Hence, the properties of these 
states are inherently sensitive to the the magnetic field 
and the details of potential profile. In contrast, the kink- 
states are of topological origin and so persist even in the 
absence of a magnetic field. Due the strong confined of 
these states, the properties of these states are remarkably 
robust and insensitive to the magnetic field (see Zarenia 
et al^foT more details). 

The LDOS of the A2 sites in the vicinity of a kink is 
shown in Fig. |11[ It was computed using states between 
an energy window between Etip and fi = O.Ot in order 
to be representative of electron tunnelling. Away from 
the interface, sharp LL plateaus easily distinguish the LL 
spectrum. The marked suppression of the LDOS on the 
other side of the interface where the parity of the bias 
is reversed is indicative of the strong layer polarizability 
discussed previously. Only when the Eup lies above the 
n — 2 LL is there any weight on the A2 sublattice in this 
region, consistent with the low energy theory. 

When Etip lies between the chemical potential and the 
n = 1 LL, there is notable enhancement of the LDOS 
close to the interface. A cut along Eup = 0.025i shown 



fests as a transfer of weight from the bulk to the inter- 
face and results in the combined downturn in the bulk 
and sudden enhancement in the LDOS close to the in- 
terface. Evidence of the strong localized states can also 
be seen when Eup lies just above the n — LL and /i — 
lies just below. Here, the tunnelling into the kink states 
is suppressed by Pauli blocking, leading to a sharp re- 
duction in the LDOS close to the interface. In contrast, 
at higher energies the reduction in the LDOS is a very 
gradual roUoff across the interface, showing no indication 
of localized interface modes. 

Evidence demonstrating the robustness of these modes 
in the presence of magnetic field suggest these modes 
may play an important role in establishing percolation 
networks responsible for the break down of the quantum 
HaU effect .1221 



C. Array of Kinks and Antikinks 

We begin by first considering the dilute limit where 
the kink and anti-kink soliton modes are well separated 
and decouple. In the the regime where £b < A, the LL 
states are well localized and respond only to their local 
environment. Hence, the dispersion relation is expected 
to consist of a series of bulk regions - identical to those 
generated by a uniform bias with alternating bias parity 
- connected by intermediate kink/anti-kink states. This 
viewpoint is confirmed by inspecting the dispersion re- 
lation of a single kink/anti-kink pair (which, neglecting 
the spurious edge states, compose the unit cell) shown in 
Fig. [12] (b). 

In the opposite regime where is > A, the dispersion re- 
lation in Fig. 12 (c) indicates the bulk LLs become disper- 



in Fig. 11 (b) shows signatures that this enhancement 
is due to the presence of strongly confined kink states. 
Namely, the confinement modes to the zero-line mani- 



sive in general. Here, the wavefunctions are spread over 
a wide region encompassing the kink/anti-kink of the po- 
tential profile, and thereby altering their properties. In 
contrast, with respect to the low energy kink/antikink 
states the magnetic field again does little more than shift 
their position in momentum space; the mode velocity is 
robust and the wavefunctions are rigid. Only at higher 
energies, close to the LL energy levels are the dispersive 
kink states modified. 

In sum, the effect of a increasing magnetic field at low 
energy is to shift the kink and anti-kink modes relative 
to each other (the anti-kink modes move right the kink 
modes move left in momentum space). Since the mode 
velocity of the both the kink and anti-kink are reversed 
in each valley, this shift pushing the band crossing points 
up in energy in one valley and down in the other. In this 
way, a magnetic field can be used to control the band 
crossing points between modes of adjacent wires. 

Further insights into the robustness of the low energy 
modes present in the continuum model of the ID inter- 
layer bias SL studied above can also be made by employ- 
ing a simple low energy theory inspired by the kink/anti- 
kink viewpoint. In the dilute limit and in the absence 
of a magnetic field each kink (anti-kink) generates an 
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FIG. 11. LDOS along the y direction (in units of d) in the 
proximity to kink in the interlayer potential bias with |AV = 
O.OStj and large B = SOT to avoid extraneous finite size of 
affects, (a) LDOS across a kink for states as function of Eup 
where /i = 0.0. (b) Cut where Eup is between the n = and 
n = 1 LL {Eup = 0.0.025t). Shift of LDOS weight towards 
the interface is a signature of the kink modes, (c) LDOS with 
Etip above n — with /i in between the n = 1 and n — 
LL {Etip = 0.042f and fi = 0.037t). Evidence of a kink state 
is given by the sharp suppression of the LDOS, as the Pauli 
blocking prevents tunnelling in to kink states. 

identical copy of the low energy modes about each val- 
ley. At each K-point in Fig. [12] (a), the two modes for 
each kink (antikink) can be linearized about zero energy 
(see Ref. |23ifor more details). An effective low energy 
hamiltonian of the SL about one of the K-points can be 
written as 

71 

+i-^np.~p;)fljp^n] (13) 

where vq is the velocity of the modes about zero momen- 
tum, and c|,^„ (/pjjji) ^'^^ operators that create electrons 
on wire n, with momentum p^ about momentum close 
to one of the two zero energy crossing points labelled p* . 
(Note, the low energy theory of the opposite valley is 
identical except that the velocity of the kink/ ant i- kink 
modes are reversed.) Representing a magnetic field by 
the Landau gauge field A = —Byx, electrons hopping 
along a given wire n see an average constant gauge field 
given by An — —BXn. In a magnetic field, the low en- 
ergy effective Hamiltonian becomes 

n 

+ (-l)"(p, + nX/el - pDflJp^n] .(14) 

These results demonstrate that the full dispersion of the 
SL consists of multiple copies a single kink-antikink pair 



FIG. 12. Dispersion of BLG in the presence of kink/anti- 
kink with \AV\ = 0.08*. (a) decoupled modes with B = OT, 
(b) decoupled modes with Asep ~ lOOd and £b = 17d, (c) 
Ai,ep = 40d and Eb = 17d, and (d) Asep = 20d and Ib = 17d. 
Note, the large spread of the wavefunctions leads to coupling 
of the zero energy modes for Asep < 40ci and the qualitative 
features remain robust for similar i's/Asep. 

dispersion repeated every A/€^ in momentum space, con- 
sistent with the numerically computation of the disper- 
sions relations provided above. Notice, that the modes in 
a given wire crosses modes in the two neighbouring wires 
symmetrically about zero energy. If the magnetic field 
is weak, the kink modes of wire n cross with the anti- 
kink modes of wire n ± 1 at an energy of ±do • If the 

modes couple, level repulsion will cause the kink/antikink 
modes of the neighbouring wire to mix and split, but be- 
cause the bandcrossing is symmetric about zero energy 
and each the each mode couples to both neighbouring 
modes equally (for a symmetric SL), a zero energy band 
remains intact. 



D. Valley Filter 

In studying a kink/anti-kink single pair in the coupled 
regime, we observe a remarkable effect not observed in a 
previous studjE3 of the kink states that used a low energy 
model that did not capture the valley degeneracy lifting. 
As the wavefunctions of the kink and antikink modes be- 
gin to overlap, a bandgap opens at the band crossing 
points. Since the size of the bandgap is determined by 
the degree of overlap between the wavefunctions, it can 
be narrowed by either increasing the kinks-antikink sepa- 
ration or by reducing the interlayer bias, which causes the 
wavefunctions to spread further from the interface. Fur- 
thermore, as the magnetic field strength increases, the 
bandgap shifts to positive energy in one valley and neg- 
ative energy in the other due to the valley asymmetry in 
the band crossing point. Remarkably, for a single cou- 
pled kink/anti-kink pair an energy window opens where 
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only unidirectional modes are present along each kink 
and anti-kink (the direction of flow can be flipped by re- 
versing the magnetic field). Moreover, all the conducting 
modes have the same valley index and so acts as a valley 
filter simil ar to those suggested along edge stated^ and 
line defect 4^^'^ when the chemical potential is tuned to 
lie in the bandgap. 

IV. SUMMARY AND DISCUSSION 

In this paper, we studied magnetic properties of MLG 
and BLG under various types of ID external potentials. 
For MLG, we identified three regimes of magnetic field 
strengths that generate distinctive features in the LL dis- 
persion and transport properties. Under a weak magnetic 
field, MLG SLs exhibit zero energy LLs whose degenera- 
cies are identical to the numbers of Dirac points present 
in the spectrum. At higher energies but still within the 
linear range of the spectrum, differences between the 
Dirac cones cause the LL degeneracy to be lifted. There- 
fore, measurements that carefully determine the degener- 
acy of the LLs can be used to probe and characterize the 
underlying anisotropic Dirac cones. We further showed 
that the diagonal conductivities show strong anisotropy, 
with conductivity in the SL direction larger than that 
in the transverse direction. Interestingly, the anisotropy 
can be reversed for an intermediate magnetic field, where 
the LLs become dispersive. This field tunable anisotropy 
may find interesting device applications. 

We then considered two types of SLs for BLG and 
again showed that analysis of the LL spectrum provides 
clear signatures of the underlying Dirac cones. However, 
the effects of a magnetic field were shown to be rather dif- 
ferent in the bilayer SL systems than on the monolayer SL 



system. Although the diagonal conductivity of both SLs 
also have strong anisotropy like in the monolayer, there 
is no field tunable anisotropy reversal. Furthermore, re- 
gardless of the type of SL there are distinct crossovers 
from nonrelativistic physics to relativistic Dirac physics 
as the SL strength is tuned. In the case of a chemical po- 
tential modulation, the zero energy LL vanishes at a crit- 
ical SL strength where the Dirac cones becomes gapped. 
As for the electric field SL, two zero energy LLs are al- 
ways present irrespective of the strength of the magnetic 
field SL strength or magnetic field. This demonstrates a 
remarkable robustness of the Dirac points inherited from 
their topological origin. 

We have also studied how a magnetic field affects BLG 
subject to a uniform bias and with bias reversing kinks 
in the potential profile. The motivation here was to gain 
a real space perspective of the LL wavefunctions in the 
quantum hall regime. Together with a low energy model, 
this study further established an intuition into the prop- 
erties of the electric field SL. Moreover, it provided valu- 
able insights into the topological kink states present in 
disorder BLG and suggests how that they may contribute 
to the breakdown of the quantum hall effect. Finally, we 
proposed that a pair of coupled kink-antikink modes sub- 
ject to a magnetic field could serve as a possible route to 
fabricate a switchable one-way valley filter. 
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